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ABSTRACT 

Scramjet  operation  in  the  lower  hypersonic  regime  between  Mach  4  and  8  is  characterized  by  what  is 
called  dual-mode  combustion.  In  this  situation  disturbances  generated  by  heat-release  in  the  combustor 
can  propagate  upstream  of  fuel  injection  to  affect  the  operation  of  the  inlet.  The  method  use  to  alleviate 
this  problem  is  installation  of  a  short  duct  between  the  inlet  and  the  combustor  known  as  an  isolator.  This 
article  describes  the  flow  phenomenon  that  exits  in  the  isolator  in  these  situations  and  presents  some  of 
the  current  methodologies  and  analysis  techniques  for  scramjet  isolator  design. 


NOMENCLATURE 


A 

area  (m2) 

V 

velocity  (m/s) 

cp 

specific  heat  (J/kgK) 

X 

axial  distance  (m) 

cf 

skin  friction  coefficient 

<t> 

equivalence  ratio 

D 

hydraulic  diameter  (m) 

T 

constant  in  mixing  curve 

fst 

stoichiometric  ratio 

Y 

ratio  of  specific  heats 

F 

stream  thrust  (N) 

0 

boundary  layer  momentum  thickness  (m) 

hpr 

heat  of  combustion  (J/kg  of  fuel) 

Pc 

combustion  efficiency 

H, 

M 

P 

total  enthalpy  (OK  basis)  (J/kg) 

Mach  number 

pressure  (Pa) 

P 

X 

density  (kg/m3) 

shear  stress  (Pa) 

Q 

heat  loss  to  the  structure  (J/kg) 

Subscript 

R 

gas  constant  (J/kgK) 

c 

core  flow 

Re 

Reynolds  number 

t 

total 

s 

T 

shock  train  length  (m) 

temperature  (K) 

w 

wall 

1.0  INTRODUCTION 

Scramjet  operation  in  the  lower  hypersonic  regime  between  Mach  4  and  8  is  characterized  by  what  is 
called  dual-mode  combustion.  Figure  1  shows  a  schematic  of  a  scramjet  powered  vehicle  operating  in  this 
way.  In  this  instance,  flow  is  compressed  by  shock  waves  in  the  forebody  and  inlet,  and  is  supplied  to  the 
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ORGANIZATION 


combustor  at  supersonic  conditions.  Combustion  of  fuel  with  the  incoming  air  generates  a  large  local 
pressure  rise  and  separation  of  the  boundary  layer  on  the  surfaces  of  the  combustor  duct.  This  separation, 
which  can  feed  upstream  of  the  point  of  fuel  injection,  acts  to  further  compress  the  core  flow  by  generating 
a  series  of  shock  waves  known  as  a  shock-train.  A  short  length  of  duct,  called  the  isolator,  is  usually  added 
to  the  scramjet  flowpath  upstream  of  the  combustor  to  contain  this  phenomenon  and  stop  it  from 
disrupting  the  operation  of  the  inlet.  In  some  engines,  the  combination  of  diffusion  in  the  isolator  and  heat 
release  in  the  combustor  decelerates  the  core  flow  to  subsonic  conditions.  In  this  instance  the  core  flow 
must  then  re-accelerate  through  Mach  1  in  what  is  known  as  a  thermal  throat. 


Propulsion-airframe  integrated  scramjet 


Figure  1:  Schematic  of  a  scramjet  operating  in  dual-mode  (NASA). 

Dual-mode  combustion  can  produce  large  pressure  levels  in  the  combustor  and  nozzle,  generating  high 
levels  of  thrust.  This  flow  is  affected  by  many  parameters,  including  the  state  of  the  boundary  layer  in  the 
isolator,  the  flow  Mach  number  exiting  the  inlet,  the  area  distribution  of  the  combustor,  and  the  position 
and  number  of  fuel  injection  stations.  As  the  separated  regions  on  the  surfaces  of  the  isolator  and 
combustor  are  seen  by  the  core  flow  as  blockage,  scramjet  engines  operating  in  dual-mode  can  be  thought 
of  involving  fluid-dynamic  variable  geometry.  At  speeds  above  Mach  8,  the  increased  kinetic  energy  of 
the  airflow  through  the  engine  means  that  the  combustion  generated  pressure  rise  is  not  strong  enough  to 
cause  boundary  layer  separation.  Flow  remains  attached  and  supersonic  throughout  in  the  instance,  and  the 
engine  operates  as  a  pure  scramjet. 

The  article  will  first  describe  the  flow  structure  that  occurs  in  the  isolator  of  a  dual-mode  scramjet.  Next,  a 
diffuser  model  applicable  to  the  analysis  of  dual-mode  scramjet  combustion  will  be  presented.  Some 
results  of  the  analysis  will  then  be  compared  to  experiments. 


2.0  FLOW  STRUCTURE  IN  A  SCRAMJET  ISOLATOR 

The  structure  of  the  supersonic  flow  in  confined  ducts  under  the  influence  of  a  strong  adverse  pressure 
gradient  is  of  interest  in  the  design  of  scramjet  isolators.  As  shown  in  the  schematic  of  Fig.  2,  the  pressure 
gradient  is  imposed  on  the  incoming  supersonic  flow  in  the  form  of  shock  waves.  If  there  were  no 
boundary  layer,  a  normal  would  form  in  a  plane.  However,  the  presence  of  an  incoming  boundary  layer 
produces  a  series  of  normal  or  oblique  shocks  that  can  spread  the  pressure  rise  over  a  length  of  many  duct 
diameters.  This  phenomenon,  known  as  a  “pseudo  shock”1  or  “shock-train”  is  characterized  by  a  region  of 
separated  flow  next  to  the  wall,  together  with  a  supersonic  core  that  experiences  a  pressure  gradient  due  to 
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the  area  restriction  of  the  separation,  forming  a  series  of  crossing  oblique  shocks  in  the  core  flow.  A 
mixing  region  also  grows  between  the  core  and  separated  flows,  balancing  the  pressure  rise  in  the  core 
against  the  shear  stress  on  the  boundary  of  the  separation.  Finally,  the  flow  reattaches  at  some  point  and 
mixes  out  to  conditions  that  match  the  imposed  back-pressure.  Being  able  to  predict  the  length  scale  of 
this  flow  structure  is  the  key  component  of  isolator  design  for  dual-mode  scramjets. 


Reattachment 


^7^77^77V^^777777777777777^^^^7777777777777777777. 

Shock  Train  '  '  Poninn 


Mixing  Region 


Figure  2:  Schematic  of  flow  structure  in  an  isolator. 


Important  work  on  this  phenomenon  was  performed  by  McLafferty2  in  the  1950’s  and  by  Waltrup  and 
Billig3  in  the  early  1970’s.  In  this  research  it  was  observed  that  for  a  given  imposed  pressure  rise  (AP/P), 
the  length  over  which  the  shock  train  spread  varied  with  DI/20' 2  and  inversely  with  (M2  -  l)(Ree)14,  where 
D  is  the  duct  diameter,  M  is  the  Mach  number  of  the  inflow,  0  is  the  momentum  thickness  of  the  boundary 
layer  and  Reeis  the  Reynolds  number  based  on  momentum  thickness.  Based  on  experiments  in  round  ducts 
with  incoming  Mach  numbers  between  1.5  and  2.7,  Waltup  and  Billig  found  that  thicker  boundary  layers 
lead  to  longer  shock  trains, and  developed  the  following  empirical  correlation  for  s,  the  distance  over 
which  the  shock  structure  (or  pressure  rise)  is  spread: 


s(M2  -l)(Reg) 

DxnGm 


1/4 


=  50 


A P 


+  170 


A P 


(1) 


Figure  3  below,  taken  from  Ref.  3,  shows  the  form  of  the  simple  quadratic  correlation  in  AP/P  in 
comparison  to  the  data  used  to  develop  the  correlation. 
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Figure  3:  Correlation  of  experimental  shock  train  data3. 


It  has  been  postulated  by  many  authors  that  the  pressure  gradient  experience  by  the  core  flow  in  the  duct 
must  be  equal  to  the  pressure  gradient  that  can  be  supported  by  shear  in  the  separated  region.  Based  on  a 
large  amount  of  experimental  data  at  different  Mach  numbers,  Reynolds  numbers  and  in  different  duct 
geometries,  Ortwerth4  determined  that  the  rate  of  pressure  rise  (diffusion)  in  a  duct  is  directly  proportional 
to  the  dynamic  pressure  of  the  incoming  flow  and  the  skin  friction  coefficient  at  the  initial  point  of 
separation  in  the  duct,  and  inversely  proportional  to  the  duct  hydraulic  diameter.  From  this  he  developed  a 
diffuser  model  for  separated  flow  in  ducts  which  can  be  expressed  as: 


dP  89  pV\ 
dx~  Dh  fo{  2  j 


(2) 


where  DH  is  the  hydraulic  diameter  of  the  duct,  Cm  is  the  friction  coefficient  at  the  initial  separation  point. 
In  essence,  this  relationship  supplies  the  ability  to  determine  a  length  scale  over  which  pressure  rise  must 
be  spread  in  a  duct.  It  will  be  used  in  the  next  section  of  this  paper  as  the  extra  equation  needed  to  perform 
quasi-one-dimensional  calculations  of  flow  properties  in  separated  ducts. 


3.0  DIFFUSER  MODEL  FOR  SCRAMJET  ISOLATOR  FLOW  ANALYSIS 

Analysis  of  the  combustion  process  in  a  scramjet  usually  involves  quasi-one-dimensional  cycle  analysis 
methods.  While  the  real  combusting  flow  in  a  scramjet  is  far  from  uniform  at  any  cross-section  throughout 
the  engine,  when  used  properly,  these  techniques  provide  an  efficient  means  of  modeling  isolator  and 
combustor  region  of  a  scramjet.  While  some  methods  simply  jump  from  the  start  to  the  end  of  the 
combusting  zone5,  the  method  presented  in  this  article  enables  prediction  of  the  pressure  distribution  in  the 
entire  region  of  the  engine  affected  by  combustion,  therefore  enabling  comparison  with  experimental 
pressure  measurements.  These  methods  follow  directly  from  the  classical  quasi-one-dimensional 
gasdynamics  presented  by  Shapiro6,  with  the  addition  of  Ortwerth’ s  diffuser  model  to  close  the  equation 
set. 

A  differential  element  of  the  separated  flow  in  a  duct  is  shown  in  Fig.  4.  Here  the  area  of  the  core  flow 
passing  through  the  duct  (Ac),  is  equal  to  the  geometric  area  of  the  duct  (A)  if  flow  is  attached,  but  is  less 
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than  A  for  separated  flow.  Fuel  and  air  are  burning  in  this  element,  and  a  friction  force  dFr  =  xwAw  is 
applied  by  the  walls,  together  with  a  heat  loss  in  the  amount  dQ.  For  simplicity  of  analysis,  the  flow  is 
assumed  to  be  that  of  a  calorically  perfect  gas  with  constant  ratio  of  specific  heats,  y,  gas  constant  R,  and 
specific  heat  at  constant  pressure,  cp.  Combustion  heat  release  is  modeled  through  the  use  of  a  heat  of 
combustion,  hpr,  and  the  change  in  total  enthalpy  of  the  flow  as  it  traverses  the  element  is: 

dH,  =  hprfstd(j)  -  dQ  (3) 

where  fst  is  the  stoichiometric  fraction  of  fuel  to  air,  and  di()  is  the  equivalence  ratio  of  fuel  that  combusts  in 
length  dx.  The  corresponding  change  in  the  total  temperature  of  the  flow  is  therefore  dTt  =  dl  lt/cp.  The 
wall  shear  stress  is  related  to  a  skin  friction  coefficient  through  xw  =  CfpV2/2,  and  Aw  =  4Adx/D,  where  D 
is  the  hydraulic  diameter  of  the  duct. 


Figure  4:  Differential  element  of  separated  flow. 


The  differential  conservation  equations  of  mass,  momentum  and  energy  for  the  element,  are  given  by: 

dp  dV  dAr 

—  + - +  — —  =  0 

p  y  a 

dp  yM1  4Cfdx  rM2  Ac  dV 2 


-  + 

P  2 


D 


2  A  V 


=  0 


(4) 


(5) 


dT  y-\  2  dV2  f 

—  +  - M — r 

T  2  V 


i+— — -m2 


dT, 


(6) 


Note  that  these  equations  all  apply  to  the  core  flow  area,  but  friction  and  heat  loss  are  based  on  the 
geometric  area.  Together  with  the  equation  of  state  for  the  gas  and  the  definition  of  Mach  number  (in 
differential  form): 


dp  dp  dT  _  q 
p  p  T 


RTO-EN-AVT-185 


10-5 


Scramjet  Isolators 


ORGANIZATION 


dM 2  dV2  dT_ 
M2  V2  +  T 


(8) 


we  have  five  equations  to  relate  the  eight  variables.  Following  Shapiro6,  area  change  (dA/A)  and  total 
temperature  change  (dTt/Tt)  are  treated  as  independent  variables.  If  the  flow  is  attached,  this  closes  the 
equation  set  and  all  1-D  axial  property  distributions  in  the  duct  can  be  calculated  as  described  by  the 
influence  coefficients  of  Shapiro.  For  separated  flows,  however,  an  extra  equation  is  needed  for  the  extra 
variable,  Ac.  This  is  supplied  by  Ortwerth’s  diffuser  relation  (eqn  2).  After  a  significant  amount  of 
algebraic  manipulation  of  equations  2,4-7,  the  following  differential  relation  for  Mach  number  is  obtained: 


d(M2) 
M 2 


( 

1  + 

V 


X 


dpi  p 
yM2  AC 
2  A 


4  C, 


+  - 


dx 

D 


A 

A 


+  - 


dT, 


(9) 


This  related  the  change  in  Mach  number  to  the  amount  of  diffusion  (dP/P),  the  amount  of  heat  release 
(dTt/Tt),  and  the  axial  distribution  of  Cf.  This  must  be  integrated  in  conjunction  with  the  following  relation 
for  Ac/A: 


d(Ac  /  A)  

dp 

f  l  +  (y-\)M2^ 

.  „  dx  (  y-\  t2\dT 

/[C  i  i  '  ]\/f-  < 

AJA 

yM2Acl  A 

P 

{  2  AJA  j 

1  D  {  2  J  Tt 

(10) 


To  determine  the  axial  distributions  of  Mach  number  and  Ac/A  in  ducts  with  specified  area  A(x)  and  prior 
knowledge  of  the  heat  release  distribution,  Tt(x)  and  Cf,  these  equations  can  be  integrated  with  a  standard 
ODE  solver  for  multiple  equations. 

This  methodology  will  now  be  applied  to  a  scramjet  vehicle  flying  at  Mach  7  and  a  dynamic  pressure  of 
50  kPa.  In  keeping  with  the  notation  of  Heiser  and  Pratt7,  station  0  is  in  the  ffeestream  flow  ahead  of  the 
vehicle,  and  a  streamtube  with  area  A0  is  captured  and  processed  by  the  engine.  Station  1  is  downstream  of 
the  vehicle  forebody  shock  and  represents  the  properties  of  the  flow  that  enters  the  inlet.  Station  2  is  at  the 
inlet  throat,  which  is  usually  the  minimum  area  of  the  flowpath,  and  the  length  between  stations  2  and  3  is 
referred  to  as  the  isolator.  Station  3  represents  the  start  of  the  combustor,  and  fuel  and  air  is  mixed  and 
burned  by  the  end  of  the  combustor  at  station  4.  The  nozzle  includes  an  internal  expansion  up  to  station  9, 
and  an  external  expansion  to  station  1 0  at  the  end  of  the  vehicle. 

Assuming  a  fairly  typical  forebody/inlet  compression,  the  properties  at  the  inlet  throat  (station  2)  will  be 
assumed  to  be  M2  =  3.60,  p2  =  50  kPa,  T2  =  650  K,  Ht2  =  2.35  MJ/kg.  The  axial  distribution  of  properties 
in  a  round  isolator/combustor  duct  with  an  initial  diameter  of  0.06m  and  a  divergence  with  area  ratio  of  2 
will  be  calculated  for  different  fuelling  levels.  The  start  of  the  isolator  is  assumed  to  be  at  x2  =  0.0  m,  and 
hydrogen  fuel  (hpr  =  120  M/kg)  is  injected  at  x3  =  0.2  m  (the  start  of  the  combustor).  The  amount  of  fuel 
that  is  allowed  to  react  with  the  air  at  a  particular  station  is  dictated  by  a  mixing  efficiency  curve, 
r|c(X),  that  takes  the  form: 


>1  =  >1 


<9X 


1  +  G9-1)X 


(11) 


where  qc  tot  is  the  combustion  efficiency  at  the  end  of  the  combustor,  X  =  (x-x3)/(x4-x3)  and  9  is  an 
empirical  constant  of  order  1  to  10  which  depends  on  the  rate  of  mixing7.  For  the  current  study  r|c  tot  was 
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set  to  0.8  at  all  times  and  a  value  of  &  =  5.0  was  used.  These  values  correspond  to  robust  combustion  in 
the  engine.  Making  use  of  eqn.  11,  the  heat  release  curve  is  therefore: 

Tt  =  Tti+ihpJsrfric ~dQ)lcP  (! 2) 

Skin  friction  was  calculated  assuming  a  constant  skin  friction  coefficient  of  Cf  =  0.002  and  heat  loss  to  the 
structure  (dQ)  was  calculated  using  Reynolds  analogy  and  an  assumed  wall  temperature  of  Tw  =  600  K. 

Given  the  limitation  of  constant  y,  R  and  cp  in  the  analysis,  eqns.  9  and  10  are  integrated  in  sections  along 
the  duct.  In  the  isolator  section  upstream  of  fuel  injection,  values  of  y  =  1.37,  R  =  287  J/kg/K  and  cp  = 
1063  J/kg/K  were  used.  In  the  combustor,  average  values  of  y  =  1.31,  R  =  297  J/kg/K  and  cp  =  1255  J/kg/K 
were  used,  so  as  to  model  the  properties  of  the  real  fuel/air/combustion  products  mixture  which  vary  along 
the  length  of  the  combustor. 

Figure  5  shows  calculated  1-D  flow  properties  in  the  duct  for  fuelling  at  an  equivalence  ratio  of  ())  =  0.5.  In 
the  isolator  section  of  the  duct  the  Mach  number  reduces  and  the  pressure  and  temperature  increase  due  to 
the  action  of  friction  on  the  duct  surfaces.  At  the  start  of  the  combustor,  flow  properties  are  recalculated  to 
be  consistent  with  the  values  of  y  and  R  used  in  the  combustor  integration,  while  conserving  fluxes  of 
mass,  momentum  and  total  enthalpy  across  the  boundary  between  the  isolator  and  combustor.  Fuel  is  also 
added,  and  combustion  along  the  duct  leads  to  a  drop  in  the  Mach  number,  an  increase  in  the  temperature, 
and  the  pressure  varies  smoothly  in  response  to  the  competing  effects  of  combustion  and  area  increase. 
The  peak  pressure  and  temperature  in  the  duct  are  P/P2  =  2.19  and  T/T2  =  2.61,  and  the  minimum  Mach 
number  is  M  =  1.48.  The  analysis  results  in  an  estimate  of  the  one-dimensional  supersonic  properties  of 
the  flow  as  it  exits  the  combustor  at  x4  =  0.5  m. 


M2  =  2.65,  p2  =  50  kPa,  T2  =  650  K,  =  1 .59  M J/kg,  4  =  0.5 


Figure  5:  Fuelling  at  <(>  =  0.5;  attached  flow  through  the  isolator/combustor  duct. 


As  the  fuelling  level  is  increased,  the  pressure  rise  in  the  duct  can  reach  the  point  where  flow  separation 
can  occur.  In  the  current  analysis  the  well  known  separation  criterion  of  Korkegi8  are  used  to  determine  if 
this  occurs.  When  this  does  happen,  an  iterative  technique  is  required  to  determine  the  axial  position  of 
separation,  after  which  the  core  flow  area  is  less  than  the  geometric  area  until  re-attachment  occurs.  It  has 
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been  found  that  in  all  instances  there  is  a  unique  position  for  the  separation  point  that  allows  the  flow  to 
re-attach  smoothly  in  the  divergent  section.  Furthermore,  if  the  core  flow  reduces  to  subsonic  conditions  in 
the  separated  region,  the  flow  re-attaches  subsonically  and  then  re-accelerates  through  a  thermal  throat  at 
an  axial  position  that  can  be  calculated  a  priori,  as  outlined  in  Shapiro6. 

Figure  6  shows  the  calculated  properties  for  a  fuelling  level  of  c|)  =  0.72.  Flere  the  pressure  rise  is  high 
enough  to  separate  the  flow,  and  the  analysis  has  been  iterated  to  determine  a  separation  point  at  x  =  0.180 
m.  The  core  flow  begins  diffusing  at  this  point  at  a  rate  dictated  by  eqn.  2,  reaching  a  minimum  area  of 
Ac/A2  =  0.822.  Combustion  of  fuel  acts  to  push  the  flow  towards  re-attachment,  which  occurs  at  x  =  0.213 
m  (downstream  of  fuel  injection)  with  M  =  2.027.  After  re-attachment  the  Mach  number  continues  to  drop 
to  a  minimum  of  M  =  1.087,  and  the  pressure  continues  to  rise  to  a  maximum  of  P/P2  =  3.51,  after  which 
the  flow  accelerates  under  the  action  of  the  increasing  area  to  leave  the  combustor  supersonically. 


x  (m) 


Figure  6:  Fuelling  at  4>  =  0.72;  separated  flow  with  a  supersonic  re-attachment. 


Figure  7  shows  the  calculated  properties  for  a  fuelling  level  of  (()  =  0.81.  Once  again  the  pressure  rise  is 
high  enough  to  separate  the  boundary  layer,  but  in  this  instance  the  increased  combustion  has  pushed  the 
separation  well  upstream  of  fuel  injection  to  x  =  0.099  m.  Furthermore,  the  Mach  number  of  the  core  flow 
has  dropped  below  Mach  1  to  form  what  is  known  as  a  “thermal  throat”.  After  separation  the  core  flow 
reduces  in  area  to  a  minimum  of  Ac/A2  =  0.822  at  the  start  of  the  combustor.  The  combustion  then  acts  to 
push  the  core  flow  to  re- attachment,  and  the  Mach  number  drops  through  M  =  1  with  re-attachment  at  x  = 
0.284  with  M  =  0.960.  The  pressure  also  peaks  at  the  point  of  re-attachment  at  P/P2  =  4.24.  The  flow  then 
re-accelerates  through  a  thermal  throat  at  x  =  0.295  m  and  leaves  the  combustor  supersonically.  There  is 
an  important  discussion  in  Shapiro6  on  the  nature  of  a  thermal  throat  and  the  manner  in  which  attached 
flow  can  pass  through  the  sonic  point. 


10-8 


RTO-EN-AVT-185 


Scramjet  Isolators 


x  (m) 

Figure  7:  Fuelling  at  <J>  =  0.81 ;  separated  flow  with  a  thermally  throated. 


A  comparison  of  Figs.  5,  6  and  7  shows  the  upstream  progress  of  the  combustion  influence  with  increasing 
heat  release.  For  <|)  =  0.81,  the  separation  is  approximately  half  way  along  the  isolator,  and  a  designer 
would  have  to  decide  if  extension  of  the  isolator  is  justified  to  allow  for  increased  combustion. 

While  it  is  recognized  that  this  analysis  involves  the  significant  assumption  of  a  perfect  gas,  it  does 
however  contain  all  the  physical  attributes  that  are  exhibited  by  real  flows.  Similar  analyses  of  combustion 
flows  using  finite  volume  techniques  and  equilibrium  chemistry  are  presented  by  Auslender  &  Smart9. 


4.0  CONCLUDING  REMARKS 

The  key  feature  of  isolator  design  is  the  choice  of  the  length  required  to  “isolate”  the  inlet  from  influences 
propagating  upstream  from  the  combustor.  Determination  of  this  length  requires  modeling  of  separated, 
diffusing  flows  in  internal  ducts.  The  diffuser  model  of  Ortwerth  has  been  implemented  here  as  part  of  a 
quasi-one-dimensional  cycle  code  for  the  calculation  of  these  flows.  Given  the  distribution  of  heat  release 
in  the  combustor  and  the  isolator/combustor  geometry,  this  code  predicts  the  length  of  the  upstream 
influence  and  hence  the  required  isolator  length. 
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